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ABSTRACT 


I 


Is 


1 


Let  X^  K  •••  K  ^(n)  ke  *he  or^er  statistics  of  a  random  sample  from  a 
population  with  density  f  and  distribution  function  F  such  that  F(0)  *=  0  . 
Let  q(t)  =  f(t)(l  -  F(t))  *  be  the  failure  rate  of  F  .  In  testing 
Hq:  q(t)  *  A  vs.  q(t)f  ,  Prosehan  and  Fyke  (Vth  Berk.  Symp.)  considered 

certain  statistics  based  on  R^,  ...,  R^  ,  the  ranks  of  the  normalized  sample 
spacings  =  (n  -  i  +  -  X^  ^),  1  <  i  <  n,  X^  *=  0  .  They  show  that 

these  statistics  are  asymptotically  normal  for  fixed  F  and  compute  the  efficacy 
of  one  of  these  statistics  for  selected  distributions.  We  show  that  asymptotic 
normality  holds  also  for  sequences  of  alternatives  approaching  Hq  as  n  -*■  » 
and  conclude  that  the  above  efficacies  yield  Pitman  efficiencies.  The  statistic 

V  ■  i  Rj  is  asymptotically  equivalent  to  the  one  considered  by  Prosehan 
and  Pyke.  If  S  *  i  log  [1  -  R^(n  ▼  1)  *]  then  the  Pitman  efficiency  of 

V  to  S  is  3/4  .  If  jfg  |  is  a  sequence  of  alternative  densities,  let 

h(t)  -  [3  log  fa(t)/30]e_o,  e±  -  (/*  h’(t)  exp(-t)dt  ]  [1  -  ifn  +  l)"1]”1,  a±  - 

-  log  {1  -  i(n  +  1)  *}  .  ’Then  £  c^D^  *9  asymptotically  most  powerful  for 
fg  )  if  the  scale  parameter  A  is  known.  However,  the  statistic  S  is 

ni 

nowhere  asymptotically  most  powerful,  although  it  is  the  asymptotically  most 

powerful  linear  rank  statistic  for  a  suitable  parametric  family  lf„  I  A 

1  ni 

comparative  study  of  these  and  other  statistics  is  given  in  terms  of  Pitman 
efficiencies,  and  Monte  Carlo  power. 
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TESTS  FOR  MONOTONE  FAILURE  RATE  BASED  ON 
NORMALIZED  SPACINGS 

by 


Pe;er  J.  Bickel 
Kjell  A.  Doksum 


1 .  Introduction  and  Summary 

Let  F  be  a  distribution  with  density  f  ,  and  let  q(t)  ■  f(t)/[l  -  F(t)] 

be  the  failure  rate  of  F  .  Tests  for  constant  versus  monotone  increasing 

failure  rate  based  on  the  ranks  of  the  normalized  spaclngs  between  the  ordered 

observations  have  been  considered  by  Proschan  and  Pyke  [10].  They  show  that 

thes'<:  statistics  are  asymptotically  normally  distributed  for  fixed  alternatives 

F  and  compute  the  ratios  of  the  efficacies  of  one  of  their  rank  tests  to  the 

best  statistics  for  Weibull  and  Gamma  alternatives. 

In  this  paper,  it  is  shown  that  asymptotic  normality  holds  also  for  sequences 

of  alternatives  JFQ  |  that  approach  the  H  distribution  1  -  exp(-Xt),  t  >  (>  , 
I  n! 

as  n  -*■  •  ;  and  that  the  above  mentioned  ratios  of  efficacies  are  in  fact 
Pitman  efficiencies. 

Let  R^,  . . . ,  Rn  be  the  ranks  of  the  normalized  spacing,  T^  ■  £  i  R^ 
and  Tj  ■  -£  i  log[l  -  R^/(n  +1)]  .  Then  T^  is  asymptotically  equivalent  to 
the  Proschan  Pyke  statistic.  It  is  shown  that  the  Pitman  efficiency  satisfies 

e(Tr  T2)  =  3/4  (1.1) 

for  all  sequences  of  alternatives  jFg  j  and  thus  T^  is  asymptotically 
inadmissible. 

Statistics  that  are  linear  in  the  normalized  spaclngs  and  asymptotically 

a 

most  powerful  for  parametric  alternatives  |F.  | ,  if  the  scale  parameter  X  is 

I  6nl 

known,  are  derived,  and  it  is  shown  that  the  rank  statistics  that  are 
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asymptotically  most  powerful  in  the  class  of  linear  rank  tests,  are  nowhere 
most  powerful  in  the  class  of  all  tests,  when  X  is  known. 

If  X  is  unknown,  studentizing  of  the  linear  normalized  spacing  tests 
which  are  asymptotically  most  powerful  for  X  known  lead  to  procedures 
which  have  only  the  same  asymptotic  power  as  the  most  powerful  linear  rank 
tests . 

Unbiasedness  is  shown  for  tests  that  are  monotone  in  the  normalized 
spacings,  and  Monte  Carlo  power  estimates  are  used  to  compare  the  various 
statistics  with  the  likelihood  ratio  tests  considered  by  Barlow  [1]. 


\ 


2 .  Tests  Monotone  in  the  Normalized  Spaejngs 

Let  X^,  . . . ,  X^  be  a  random  sample  from  a  population  with  a  continuous 
distribution  F  satisfying  F(0)  =  0  ,  and  let  0  ■  ^(o)  <  ^(1)  <  •••  <  ^(n) 
be  the  order  statistics.  The  normalized  sample  spacings  D^,  ...»  are 

defined  by  Di  =  (n  -  i  +  1)  (X^  -  X^^),  i  ■  1,  ....  n  ;  and  denotes 

the  rank  of  among  D^,  .  The  problem  is  to  test 

H  :  "-log[l  -  F(x)]  =  Xx  on  [0,  °)  for  some  positive  constant  A"  against 
Hj  :  "-log(l  -  F(x)]  is  convex  and  not  of  the  form  Ax  on  [0,  00 )  Note 
that  the  only  distributions  satisfying  F  are 


K.  (x)  «=  1  -  e  ;  x  >  0,  A  >  0 


(2.1) 


Under  H  it  is  well  known  that  (D,  ,  ....  D  )  has  the  same  distribution  as 
o  in 

a  random  sample  from  a  population  with  distribution  ,  while  under  the  alter¬ 
native  there  is  a  downward  trend  In  the  sense  that  P(D^  >  D^)  <  whenever 

i  >  i  (see  [10]).  One  thus  defines  a  test  $  ■  $(D. ,  ...»  D  )  to  be  monotone 

in 

in  the  D's  if 


♦  (D{,  ....  D')  <  ♦(D. . D  )  for 

i  n  -  i  n 

all  (D, ,  ...,  D  )  and  (Di  ,  D')  such  (2.2) 

in  in 

that  i  <  j  and  implies  . 

Following  van  Zwet  [9],  one  defines  a  distribution  to  have  a  more  slowly 

increasing  failure  rate  than  F  ,  written  F^  >  F  ,  if  F^*  F  is  convex.  Here 
f"1  F  is  defined  by  P^'1  F(X)  <  x  |  f)  «  F^x),  x  >  0  . 

Theorem  2.1: 

Monotone  tests  have  monotone  power ,  i.e.’,  if  $  is  a  monotone  test  and  if 


i 


A 


4 

| 

i 

r 
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F.  >  F  ,  then 
1  c 

E($  |  F1)  <  E($  |  F)  .  (2.3) 

Proof : 

Since  F^1  F  is  increasing,  «=  F^1  F(X^)  is  the  ith  order 

statistic  in  a  random  sample  from  a  population  with  distribution  F^  .  Let 

=  (n  -  i  +  1)  (X'^  -  X'^^),  i  =  1 . n  .  Since  F^  F  is  convex, 

i  <  j  and  l)J  >  Dj  implies  .  From  (2.2)  one  obtains 

4>(Dl  ,  ....  D')  <  *(0.,  ....  D  )  ,  (2. A) 

i  n  -  l  n 

and  (2.3)  follows  upon  taking  expectations  in  (2. A). 

Note  that  a  test  <{>  is  similar  if  E(<f>  |  K^)  is  independent  of  A  . 

Thus  all  rank  tests  are  similar. 

Corollary  2.1: 

All  monotom  tests  are  similar  and  unbiased. 

Proof ; 

Similarity  follows  since  (2.2)  implies  $(D, . D  )  =  <J>  (AD- ,  ....  AD  )  . 

ini  n 

Unbiasedness  follows  by  letting  F^  *  in  Theorem  2.1. 

Corollary  2.2: 

If  -c^U)  ond  J^d)  are  nondecreasing  in  i  *  1,  ...,  n  ,  then  the 
test  that  rejects  when 


n 


l 

i=l 


Cn(1)  W  2  C1 


I 
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has  monotone  power  and  ie  unbiased. 


Proof : 


Using  the  notation  of  Theorem  2.1,  define  R^  to  be  the  rank  of 
among  D^,  .  Then  i  <  j  and  >  R^  implies  R^  >  R^  .  From 

Corollary  2  of  Lehmann  [9]  one  obtains 


n  n 

l  c  (i)  J  (Rl)  <  l  c  (i)  J  (R  )  . 
n  n  i  -  n  n  i 


(2.6) 


It  follows  that  the  test  is  monotone,  and  Theorem  2.1  and  Corollary  2.1  apply. 
Remark  2.1: 

The  results  of  this  section  can  be  used  to  obtain  bounds  on  the  power  of 
monotone  tests.  For  instance,  if  8(V^,  F2)  denotes  the  power  of  the  Proschan 
Pyke  statistic  =  "number  of  pairs  (i,  j)  with  i  <  .1  and  >  D^" 
for  the  Weibull  distribution  F_  ,  then  the  power  satisfies  B(V  ,  F)  >  8(V  ,  F0) 
for  all  distributions  F  such  that  F2^F  is  convex.  Fj)  *8  8*ven  in  [!)• 

Remark  2.2: 

Barlow  and  Proschan  [13]  have  shown  that  if  c  (1)  >  ...  >  c  in)  .  the 

n  -  -  n  * 


test  that  rejects  when 


I  D 

i-=l  n  1 
n 

l  h 

i*=l  1 


is  similar  and  unbiased. 


r 
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3 .  Asymptotically  Most  Powerful  Tests 

Let  {f  .  :  0  >  0,  A  >  0}  be  a  class  of  densities  such  lhat 
o ,  A  — 

f  ,  (x)  *=  A  exp  (-Ax),  x  >  0  ,  and  such  that 

O,  A  - 

hA(x)  “  a!  108  fe,A(x)  I  e  =  o  (3,1) 

exists. 

We  suppose  A  is  a  scale  parameter,  i.e., 


Vx<>x  i  '>  ■  p8,iw  : tJ  • 


(3.2) 


For  test  "0  =  0"  versus  ”0  >  0"  ,  the  locally  Most  powerful  test  rejects 
"0  »  0”  for  large  values  of 


Tn  "  •>>  I  \<«1>  • 

We  will  consider  sequences  of  alternatives  If.  ,1  such  that 

‘  VM 


(3.3) 


lim  n  0  ■  b  for  some  0  <  b  <  <*>  . 

n 


(3.4) 


A  sequence  If  .1  is  said  to  be  contiguous  to  f  (in  the  sense  of  LeCam  - 
I  n*  I  °* 

Hajek)  if  for  any  sequence  of  random  variables  R^X^,  ...»  X^) ,  R^  -*  0  in  Po 

probability  implies  Rn  0  in  P0  probability,  where  P0  denotes  the 

n 

probability  distribution  of  X.,  ...,  X  if  f  is  true.  The  following 

A  II  W|A 

condition  implying  contiguity  for  sequences  as  in  (3.4)  can  be  obtained  from 
LeCam  [ 8 ] . 


(a)  3f0  ^(x)/30  i  0  whenever  f0  ^(x)  >  0  .  (3.5) 

9  9  Oo 

(b)  For  some  6  >  0  and  all  0  c  [0,  6],  0  H(0)  *=  jf  [ 3f _  ,(x)/30]^ 

_1  0  6’X 
[fn  ,(X)]  dx  <  ®  t  and  H(0)  is  continuous  in  0  . 

0,  A 
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and  H(6)  is  continuous  in  6 


(c)  ii«|h|^o  J jh  [fe+htx(x)  "  fe,x(x)]  "  *‘Pfe,A<x)/3e]  f8^x(x)j  dx  “ 


-  0  for  8  c  (0,  6)  . 


An  easy  sufficient  condition  implying  (3.5)  (b)  and  (c)  is 


/•“f  l[  -<M-< -]  “e.A 


(x))  *:0<8<6|dx<« 


(3.6) 


for  some  A  >  0  . 

It  also  follows  from  LeCam's  work  that  under  condition  (3.5)  we  have 

-  I1”*  -  lo«  -  ° 

in  P  and  hence  in  P.  probability. 

O  w  __ 

n 


(3.7) 


According  to  Wald  [12],  a  sequence  of  level  a  tests  (8  )  is  said  to  be 

n 

asymptotically  most  powerful  (A  known)  if 

U.n  .op  (£,_,(♦„)  -  *,,»<♦„>  > •  >  «.  »,.»<*„>  ;  «>  - « • 

Using  contiguity  we  can  prove  Wald's  [12]  main  theorem  under  weaker  conditions. 


Theorem  3.1: 


Suppose  ( 3.S )  holds  and , 


o'1*  l  h.(X  ) 

1  A  1 


(3.9) 


J* 


in  Pfi  .  probability  if  n  8  ■*  •  .  Then  the  sequence  of  level  o  teste 

v  ( A  fl 


8 


n 

i  which  reject  for  suitably  large  values  of  \  h  (X  )  is  asxpnptot ically 

n  i-i  A  1 

most  powerful. 

Proof: 

The  result  is  an  easy  consequence  of  the  Neyman  -  rearson  Lemma,  (3.7) 

and  (3.9)  .  || 

It  is  easy  to  see  that  the  are  locally  most  powerful. 

Our  next  result  employs  Theorem  4.1  ot  the  following  section.  We  need 

the  notation 

•0 

o*(hj)  -  /  h*(x)  e”x  dx  ,  and  (3.10) 

Sn  -  n  **  l  »(^j)  (Dj^  “  1  *)  ,  where  (3.11) 

a(u)  -  (1  -  u)  1  /  h!(x)  e~x  dx,  0  <  u  <  1  .  (3.12) 

-log(l-u)  1 

Theorem  3.2: 

If  (3.4)  and  (3.5)  hold t  and  if  h^  of  (3.1)  satisfies  the  conditions  of 
Theorem  4.1 ,  then  XSn  o  *(h^)  has  under  Hq  asymptotical ly  a  standard  normal 
distribution.  Moreover ,  the  test  that  rejects  for  large  values  of  is 

asymptotically  most  powerful. 

Proof: 

Since  X  is  scale  parameter,  and  since  Sq  is  linear,  one  obtains 

P„  , (XS  <  x)  ■  P.  , (S  <  x)  .  Thus  one  may  assume  without  loss  of  generality 
6, A  n  -  s,i  n  — 

2 

that  X  ■  1  .  Since  o  (hj)  is  the  asymptotic  variance  of  T^  ,  the 
asymptotic  normality  follows  at  once  from  Theorem  4.1.  By  Theorem  4.1,  the 
test  rejecting  for  large  values  of  T^  is  asymptotically  most  powerful.  Thus 


9 


it  is  enough  to  show  that  T  -  >  0  in  probability.  However,  this 

n 

follows  at  once  from  the  contiguity  condition  (3.5)  and  Theorem  4.1.  || 


Remark: 

Theorem  3.2  and  the  contiguity  condition  implies  the  asymptotic  normality 

of  S  also  under  sequences  of  alternatives  for  which  (3.4)  is  satisfied,  i.e., 
n 

if 


GU 

Pn^V  =  f  h!^x)  f0  ^  M  dx  > 


(3.13) 


then  [AS  -  u  (h,)]/o(ht)  has  asymptotically  a  standard  normal  distribution, 
n  n  1  1 
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A  .  Linear  Approximations  to  Locally  and  Asyupt  o  t  ica  1 1  y  most  Powerful 
Statistics 

As  we  have  seen,  locally  and  asymptotically  most  powerful  test  statistics 
for  parametric  alternatives  are  of  the  form 


i=l 


h(X^)  =  n 


-H 


i=l  u; 


for  some  function  h  on  [0,  .  In  order  to  compute  asymptotic  of f icienoio 

we  will  need  the  following  result  which  gives  conditions  under  which  can 

be  approximated  by  a  statistic  linear  in  the  D^'s  »  and  given  by, 


S 

n 


(A.  2) 


where, 


a(u)  *=  (1  -  u)  1  J  h'(x)  e  X  dx,  0  <  u  <  1  . 
-log(l-u) 


Theorem  A.l: 


Let  h  be  any  function  such  that 


(i)  h'  is  continuous  on  (0,  °°)  , 


oo 

/ 


(ii)  /  h(t)  e  d't  =  0,  0  <  /  h^(t)  e  C  dt  <  , 


OO 

/ 


0  0 

(iii)  either  one  of  the  following  holds 


(a)  h'(-log(l  -  u)),  0  <  u  <  1  ,  satisfies  assumption  E 

of  Chernoff,  Gastwivth  and  Joh  .  (1967 )  and 


m 

/ 


h’(v) 


v 

2 


(1  -  e  v)  dv  <  •  . 


(b)  h'(t)  changes  sign  only  a  finite  number  of  times  as 

t  -*  0  or  •  and  h'  does  not  vanish  infinitely  often. 

Moreover,  suppose  that  the  X's  have  the  exponential  density  exp(-x),  x  >  0 
Then  Tn  -  Sn  tends  to  zero  in  probability. 

The  proof  is  deferred  to  the  appendix;  Section  7. 


Remark: 

It  may  also  be  shown  that  If  (i),  (ii)  and  (Hi)  (b)  hold  and  X  ■  1, 

then  E-  .(T  -  S  )2  -  0  . 

u,  l  n  n 


Let  Pitman  asymptotic  efficiency  be  defined  as  usual  (e.g.,  Hodges  and 


Lehmann  [6]).  In  this  section  we  shall  show  the  asymptotic  normality  of 
statistics  of  the  form 


W  =■ 
n 


-] 

c.  -  n 


n 


(5.1) 


and  compute  their  Pitman  efficiencies  with  respect  to  the  asymptotically  most 
powerful  statistics  of  Section  A.  Let  h^  be  as  defined  in  (3.1)  and  let 

Oo 

ai  "  a(n+l)  where  a(u)  *  (1  _  u)  1  J"  h^(x)  e_X  dx  .  Then  (5.2) 

-log(l-u) 


S 

n 


i 


x'1) 


(5.3) 


is  the  asymptotically  most  powerful  statistic  of  Section  3. 

Theorem  5.1: 


If  (S.4)  and  (3.5)  hold  ,  if 
of  this  paper,  and  if  { c  ^ }  and 


h^  satisfies  the  condition 
I J (7^1) |  satisfy  the  cond 


s  of  Theorem  4.  1 
i lions  of  Theorem  4. 


1 


of  hajek  [5],  then  for  the  sequence  of  alter>iatives  given  i?i  (3.4) 

[W^  -  y(Wn)]/o(Wn)  has  asymptotically  a  standard  normal  distribution,  where 


ai(ci 


[-log(l  -  u)  -  l]du  ,  and 


(5. A) 
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Then  Theorem  5.1  yields 


a 


.-vr  *w**j  » 


Corollary  5.1: 


Under  the  conditions  of  Theorem  5.  J ,  the  Pitman  asymptotic  relative  efficiency 


of  W  to  S  it 
J  n  n 


e(W,  S)  =  Cor  (J(U),  -log(l  - 


r  ?  Vn(a) 

U))  lim  Cor  (a,  c)  - 

n  ,  n  „ 

n-*°  It-  2 


I  V  2 

n  J  ai 


(5.10) 


where  U  is  an  unifom  (U(0,  1))  random  variable. 

It  is  clear  from  (5.9)  that  one  obtains  the  most  efficient  linear  rank 

statistic  by  taking  J(u)  ■=  -log(l  -  u)  and  c^  =  a^  .  Note  that  the 

choice  of  J  is  independent  of  the  alternative  densities  {f  }  .  Moreover, 

b  » A 

e(W,  S)  is  less  than  one  when 


lim(a) 2 


00 

/x"i 


(x)  e  X  dx  «=  Cor2(X1,  hi(X^))  t  0  . 


(5.11) 


This  is  equivalent  to  having  X^  =  n  £  X^  correlated  with  the  locally  most 

powerful  statistic  .  It  will  be  shown  in  the  proof  of  the  next  result  that 

Cor(X  ,  T  )  <  0  when  the  failure  rate  is  increasing.  Let 
n  n  ° 


q  (x,  6)  =  f  (x)  [1  -  Ffl  (x)]  ,  x  >  0  , 


(5.12) 


denote  the  failure  rate  of  F, 


We  assume  in  what  follows  that  3q^(x,  0 ) / 9 0 


exists  and  is  continuous  in  (x,  9)  for  0  <  9  <  5  .  x  >  0  ,  some  6  >  0  .  Let 


L(x)  *=  3q  (x,9)/30 


Increasing  failure  rate  clearly  implies 


6  *=  0 


(5.12) 


1 
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L(x)  >  0  (5.14) 

for  all  x  >  0  .  We  have 

Theorem  5.2: 

Suppose  the  conditions  of  Theorem  5.1  hold ,  that  (5.14)  is  satisfied  with 
strict  inequality  for  x  in  some  set  of  positive  measure  and  that t 


L(t)  dt  <  «  . 


(5.15) 


Then  each  linear  rank  statistic  of  the  form  (5.  1)  is  inefficient ,  i.e.t 


e(W,  S)  <  1  . 


(5.16) 


Proof : 


Without  loss  of  generality  let  X  -  1  .  Now, 


f e ^ i <x)  ■  q1(x,e)  exp  | 0>  dt 


(5.17) 


Since  3q^(x,  0)/30  is  continuous  in  (x,0)  we  have  differentiating  under 
the  Integral  sign. 


h^(x)  -  L(x)  -  /  L(t)  dt  . 


x 

i 


(5.18) 


Using  Fubini's  theorem  and  (5.15)  we  get  if  X  *  1  , 


cov 


(Xr  h^(Xj^))  mf%  e  *  dx  ■  -  £ i e  *L(x)  dx  <  0 

0 


(5.19) 
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by  our  assumptions.  | 


Example  5.1; 


An  alternative  for  which  there  Is  an  efficient  rank  test  is  provided  by 

f  . (x)  -  (1  -  0)"1  [1  +  0(-2x  +  4  x2)J  e"X,  6  <  4  .  (5.20) 

8,1 

For  this  density,  the  failure  rate  is  not  monotone.  It  is  easy  to  see  that 
Cor(X^,  h(X^))“  0  .  The  efficient  rank  test  rejects  for  large  values  of 


|  [lo*  (‘  -7k)  +  1]  [-*•«(*  -  ^r)] 


(5.21) 


We  conjecture  that  under  the  hypotheses  of  our  theorem  there  is  no  sequence 
of  rank  tests  which  is  asymptotically  most  powerful  in  the  sense  of  Wald.  Our 
best  approximation  to  this  result  is  Theorem  5.3.* 


Definition; 


We  say  a  sequence  of  statistics  M  (X. ,  . . . ,  X  )  is  efficient  if 

n  l  n 

...  cor.  (Cn.  £  jl=g  -  log  f0i>Cy)j  -  1  (5.22) 


for  every  sequence  of  alternatives  0r  satisfying  (3.4).  By  as.  con  we  refer 
to  the  asymptotic  correlation  computed  under  Ho  . 

It  is  easy  to  see  that  if  (4.15)  holds  snd 


R  -►  « 
n 


(5.23) 


M  . 

in  P  probability  if  n  0  ■+  •  .  then  an  efficient  sequence  of  M  s  esn 

6  o  *  n 

n 

used  to  construct  asymptotically  most  powerful  tests  for  every  level  a  . 


be 


E.  Torgersen,  using  general  considerations,  has  proved  that  our  conjecture  is 
implied  by  Theorem  5.3.  A  more  direct  proof  using  the  special  structure  of  this 
problem  may  also  be  given. 
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Theorem  5.3: 


Under  the  assumption  of  Theorem  5.2  no  efficient  sequence  of  rank  statistics 


exists. 


We  need  an  elementary  lemma. 


l.emma  5 .  A : 


Let  be  a  sequence  of  random  variables  such  that  U  ->  0  . 

.'k') 

Denote  by  IT  J  the  random  variable  equal  to  if  (u^J  <  k  and  equal  to 

k  sgn  Un  otherwise.  Then  there  exists  ,  such  that , 


ri 


+  o  . 


(5.2 A) 


Proof : 


Let  kfl  be  the  largest  k  such  that  <  n  * 

Clearly  k^  t  °°  .  |  | 


Proof  of  Theorem: 

Suppose  Mn(R^,  ....  R^)  is  an  efficient  sequence  of  rank  statistics.  Let 

-Jj 

0n  ■  n  ,  X  ■  1  .  Without  loss  of  generality  suppose, 

M„  "  L  [l08  fV1<V  '  108  ,o.i°ti)]  ' 15  Eo,i(h?(xi>)  *  0  (5-25) 


in  P  probability.  By  (3.7)  and  (3.12) 
n 


M  -  S  -*>  0 
n  n 


In  Pq  probability.  By  the  lemma  there  exist  k^  such  that 


(5.26) 


T  l<^n> 

uh-vj 


(5.27) 


1 


From  (5.27)  -  (5.29)  we  conclude  that 


Eo,i(sn  -  Eo,i<s„  I  Rr  ••••  V>  •  0 


Un  ‘  E0,l<sn  I  R1 . V  ‘  "  'I  a(n+]  )  K0.  l(X(K1>) 


is  a  statistic  of  the  form  (5.1).  Theorem  5.2  and  (5.30)  are  therefore 
incompatible.  | | 

We  conjecture  that  rank  tests  which  reject  for  large  values  of 
Eq  2 I  •••»  Rn)  are  asymptotically  most  powerful  among  all  rank  test 

for  the  given  model  whatever  be  X  . 

It  Is  not  difficult  to  see  by  using  the  method  of  Hoeffding  [7]  that  th 
locally  most  powerful  rank  test  rejects  the  null  hypothesis  for  large  values 
the  statistic 


Wn  "  E0, 1  <Tn  I  R1 . Rn) 


"  jiEoa,,jyiX(«j)<"-J  +  1>7J- 


By  the  remark  following  Theorem  4.1  since 


E 


Ur  and  are  asymptotically  equivalent.  This  lends  further  support  to  our 

conjecture. 

These  rank  statistics  are  of  course  usable  even  if  X  is  unknown,  the 
situation  which  primarily  concerns  us,  while  the  optimal  statistics  of  Section  A 
depend  on  X  and  do  not  lead  to  similar  tests.  If  we  use  the  method  of 
Barlow  [1]  and  Nadler  and  Eilbott  [1A]  and  consider  studentized  statistics  of 
the  form 


S* 

n 


(5.3A) 


it  is  not  difficult  to  show  that  for  any  fixed  X  under  the  null  hypothesis, 


S* 

n 


-4 


*  jj  (a(^x)  -  ;K  - 


x"1) 


(5.35) 


in  probability  and  hence  by  (5.6)  the  best  studentized  test  of  this  form  is 

asymptotically  equivalent  under  contiguous  alternatives  to  the  rank  test  based 

on  M  (R, ,  ....  R  )  where, 
n  l  n 


W 


Rn> 


n 

-  y 

i=l 


■(£) 


log 


(-5r). 


(5.36) 


We  conjecture  that  the  asymptotically  most  powerful  linear  rank  tests,  asymptotically 
equivalent  to  the  studentized  asymptotically  most  powerful  linear  spaclngs  tests, 
are  in  fact  asymptotically  most  powerful  among  all  level  a  tests  which  are 
similar  for  the  hypothesis  Hq  ,  X  unknown. 


Finally  we  remark  that  the  rank  statistic  whose  Monte  Carlo  power  is  studied 


20 


in  Section  6  is  not  Un  but  an  asymptotically  equivalent  simplex-  form  given  by, 


n 


suggested  by  the  discussion  following  Corollary  5.1. 
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6.  Applications  of  the  Theory  and  Monte  Carlo  Results 

The  results  of  the  previous  sections  will  now  be  applied  to  specific 
alternatives  and  specific  statistics.  The  weights  c^  of  (5.1)  11  be  of 

the  form 

Ci  "  C(n+l)  ^°r  some  funct*on  c  on  (0,  1)  (6.1) 

and  the  efficiency  (5.10)  will  be 

e (W,  S)  -  Cor2(J(U),  -log(l  -  U))  Cor2(a(U),  c(U))  Var^ W3  (6.2) 

E(a  tU)) 

where  U  is  an  uniform  (U(0,  1))  random  variable. 

The  statistics  to  be  considered  are: 
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s3  -  f  -Sf-M1  -  ir)]i  "i 
h  "  I  sfeh')!;10^  ■  si)] 

oo 

where  g(t)  =  (1  -  t)  *  J  x  3  e  X  dx  . 

-log(l-t) 

Large  values  are  significant. 

The  alternative  densities  to  be  considered  are  listed  below  for  X  *  1  ; 
to  obtain  the  general  form,  make  the  transformation  f  (x)  -+  Xf  (Xx)  . 

f *1J(x)  *=  [1  +  0(1  -  e  x)]  exp  {-[x  +  0  (x  +  e  x  -  1)]}  (Makcham) 

f^2)(x)  *=  (1  +  0x)  exp{-(x  +  *s0x2)}  (linear  F.R.) 

f<3)(x)  =  (1  +  0)  x6  exp{-x{1+0))  (Weibull) 

ff^(x)  «  (x0  e  X]/r(l  +  6)  (Gamma) 

for  each  density,  x  >  0,  0  >  0  ;  and  the  null  hypothesis  is  obtained  for 
0=0.  Each  of  these  densities  have  increasing  failure  rates  for  0  >  0  . 

Note  that  ff2^  has  the  linear  F.R.  (failure  rate)  1  +  0x  ,  while  ff3^  has 
the  failure  rate  1  +  0(1  -  e  X)  . 

From  Theorem  3.1  it  follows  that  is  asymptotically  most  powerful  for 

ff*^  when  X  is  known,  i  ■  1,  2,  3,  4  .  Theorem  5.1  implies  that  W  is 
asymptotically  most  powerful  for  in  the  class  of  linear  rank  statistics, 

i  *  1,  2,  3,  4 . 

Wq  is  asymptotically  equivalent  to  the  Proschan-Pyke  statistic  (see 

Remark  3.1).  It  is  uniformly  improved  asymptotically  by  .  We  now  list  the 
efficiencies  of  to  the  asymptotically  most  powerful  statistic  of  Section  3. 


23 


The  efficiencies  are  given  in  general  as  functions  of  h(x)  > 

[3  log  f  (x)/30)  .  .  They  are  always  independent  of  X  . 

0,1  0CU 


hjCx) 


*<y 


J  h'(x)  e  X(!5X  -  1  +  e  X)dx 

.0  _ . 


J  h2(x)  e  Xdx 


m  CO 

ll 


9  I  I  htx)  e  X(%x  -  3/2  +  2e  X)dx 


n2 


J  h2(x)  e  Xdx 


e(W1)  =  -  e(WQ) 


e(W2) 


[?>• 


(x)  e  X(%x2  -  x)dx 


I 


/  h  (x) 
0 


-x , 
e  dx 


e(W3)  - 


/  h^(x)  e  X  x(log  x  +  y  -  l)dx 

.0 


-|2 


-  2  Y  +  l|  /  h2(x)  e  Xdx 


«<V 


i  ~x 


(1  -  t)  1  (a(t)  -  a)dt  dx 


(■£  -  l)  /  h2(x)  e  xdx 


where  y  ■  .5772  is  Euler's  constant, 


,-l 


-t 


a(t)  ■  (1  -  t)  f  h'(t)  e  fcdt  ,  and  a  •  J  a(t)dt  . 

-log(l-t) 


l 


As  remarked  in  Section  5,  if 


I 
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S*  -  Sj/n'1  ^  D}  -  n*5  ^  a(^).  i  *  1 


*  1,  2,  3,  4 


(6.3) 


are  the  studentized  linear  spacings  statistics,  then  S*  have  the  same  effi  '^cy 
as  ,  i.e.,  e(S*)  »  e(W^)  i  «=  1,  2,  3,  4  .  These  efficiencies  are  given  in 

Table  6.1. 

Tables  6.2  deal  with  a  Monte  Carlo  study  of  the  powers  of  the  and 

statistics  for  the  linear  F.R.,  Weibull  and  Gamma  distributions. 

From  the  tables,  the  following  conclusions  are  apparent:  1)  The  rank  tests 
are  uniformly  less  powerful  than  the  corresponding  studentized  linear  spacings 
tests.  2)  Of  all  the  tests  considered  in  the  above  tables,  the  total  time  on  test 
statistic  S*  is  generally  best  both  on  the  basis  of  asymptotic  efficiency  and 
Monte  Carlo  power.  3)  Barlow  (1]  has  shown  that  for  Weibull  and  Gamma  alternatives, 
S*  is  much  better  than  the  1FR  likelihood  ratio  test  and  slightly  worse  than  the 
IFRA  likelihood  ratio  test. 

Table  6.2  also  gives  the  Monte  Carlo  power  of  the  test  which  rejects  for 
large  values  of 

n 

l  (1  -  x  )  log  x 
i-1  1  1 


I 


Using  the  results  of  Section  3,  we  see  that  this  test  is  the  studentized  version 
of  the  locally  most  powerful  test  for  the  Weibull  distribution.  The  tabic  shows 
the  nonrobustness  of  this  test. 

The  last  row  in  Table  6.2  gives  the  factor  by  which  the  efficiencies  have 
to  be  multiplied  in  order  to  obtain  the  efficiencies  with  respect  to  the  best 
linear  rank  tests  (or  the  best  studentized  linear  spacings  tests). 


Rank  efficiency  factor  4  2  - -  -  1.2235  — 5 —  -  2.5505 


TABLE  6.2:  Monte  Carlo  power  based  2,000  trials  for  samples 
distributions  for  significance  level  a  m  .01,  . 


J 
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Table  6.2  shows  that  for  n  =  10  and  the  distributions  considered,  the 
rank  statistics  are  consistently  less  powerful  than  the  corresponding 

studentized  linear  spacings  statistics  S*  .  This  led  us  to  compute  Table  6.3 
in  which  n  =  30  and  the  Weibull  distribution  is  considered.  This  table  shows 
that  for  small  6  ,  there  is  no  difference  in  the  powers  between  the  total  time 
on  test  statistic  S*  and  its  rank  counter  part  ,  while  for  larger  0  ,  the 
total  time  on  test  statistic  is  again  better.  However,  as  expected  from  the 
asymptotic  theory,  the  difference  in  power  seems  to  be  decreasing. 


TABLE  6.3:  Monte  Carlo  power  based  on  1,000  trials  for  samples  of  size  30  from  the  Veibull  distribution,  for 
significance  levels  a  -  .01,  .05  and  .10  . 
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7.  Append  ix 


Wo  use  the  notation  of  Section  4.  Heuristical ly  our  argument  is  very  simple 
and  is  essentially  the  same  as  that  used  in  [3]. 

_u(  n  ) 

_  *1  \  It.  /V  \  L’f  I.  /U  \l( 


Tn  =  "  j  lmJ  lh(X(i))  -  Klh(X(i))1j 


;  l  h(X(i))  "  h(E(X,n) 


U  =  i 


(i) 


•  V" 


n  (J,  h,(K(X(i)))  lX(i)  "  E(X(i))]! 


(7.1) 


i 

-  l  h’(-log(l  -  -M)  l  (I)  -  l)(n  -  j  +  l)-1 

i-1  j=l  J 


under  .  From  the  last  approximate  identity  our  result  follows.  The 

justification  of  these  approximations  poses  some  minor  technical  difficulties. 

We  proceed  with  the  proof  of  the  theorem.  We  show  that  i,  ii  and  iii  b)  suffice. 
Let , 


J^t)  *=1  6  <  *-  <  1  —  6 

*  0  otherwise 


(7.2) 


From  Theorem  3  of  [3]  it  follows  that,  if  T^  =  n  5  )  J  (—77)  h(X...)  then, 

n  ^  y  6\n+l/  U1 


n 


T6  =  n^  p{  +  n'15  l  a  (J  (n  +  1)_1)(D,  -  1)  +  o  (1) 
n  n  j  =  l  J  P 


(7.3) 


where 


1 

7 


a6(s)  *=  (1  -  s)  j  J6(t)  h'(-log(l  -  t) )dt 


(7.4) 
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6 

Un  ^ 


1=1 


J6^‘)  h«-1°C<1 


i(n  •+  l)"1)) 


(7.5) 


and  o  (1)  as  usual  denotes  a  remainder  converging  to  0  in  prob.nl)  i  1  i  ty . 

P 

Now  let 


and 


nl 


l 

j  <  (n+1) 6 


{ h (X (  } )  - 


K(h(X(j))} 


R 


-4 


n2 


l 

j  > (n+1) (1-6) 


(h(X(j)) 


r:(h(x(l)))l 


Then, 


(7.6) 


(7.7) 


(7.K) 


T 


n 


T6  +  R{. 
n  ill 


(7.9) 


We  begin  with. 


Lemma  7.1: 

If  6  -*  0  as  n  -*■  00 

n 


in  probability. 
Proof : 


(7.10) 


(7.11) 
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a  (s)  |2  <  (1  -  s)'2  f |  h'(-log(l  -  t))  |  dt  for  6  >  0  .  (7.12) 


Moreover , 


(1  -  s)  I  /  I  h ' (-log  (1  -  t ) )  |  d t  1  ds  <  ®  . 


To  see  (7.13)  note  that  the  left-hand  side  equals 


(7.13) 


r,.  »-2  / 


a .  s,-2/  |  h'(-log(l  -  t))|dt  £  |  h'(-log(l  -  v))  |  dv  ds 


(7.14) 


oo  n 

■//■  *■ 


(u)  ||  h' (r)  |  e  U(1  -  e  r)  dr  du 


after  some  standard  arguments.  Now  by  c)  there  exists  5  >  0  such  that 

h'(x)  has  constant  sign  for  x  <  6,  x  >  6  *  and  sup  .  |  h'(x)  |  <  •  . 

6  <  x  <  6 

Then 


6"1  u 


(u)  | |  h' (r)  |  e  U(1  -  e  r)  dr  du  <  ®  . 


6  6 


Suppose  for  simplicity  h'(x)  >0  x  <  6,  h'(x)  <0  x  >  6 
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6~i  6 


/  /, ,,, 

6  0 

y1 

K  /  |  h' (u) 


h'(r)  |  e  U(1  -  e  r)  du  dr 


6 

“/h' 


c  f  h(r)(l  -  e  r)  dr  du 
"0 


1 


h’(u) 


(7.15) 


h(r)  c  r  dr  +  h(6)(l  -  c  ^))  du  <  ® 


Continuing, 


since  /  e  r  |  h(r)  |  dr  <  00  . 

Jo 


6  u 

//“• 

0  0 


(u)  ||  h'(r)  |  e  U(1  -  e  r)  du  dr 


Lm  2  /  h' 
►0  / 


lim  2  J  h'(u)  e  U  (h(u)  -  h ( A ) }  du 
X- 


-f 


U 

/ 


c  Vh ' (v)  dv 


du 


du 


(7.16) 


/•- 


lim  |  e  dh^(u)  -  2h(A)  /  h'(u)  e  U  du 


■/- 


h'(u)  e  du 


6 

•/ 


h^(u)  e  U  du  -  h^(6)  e  ^ 


5 

/•" 


h(u)  du 


h( 5)  e 


-6 


2 
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+  lim  - 
A  >0 


h 2 ( A ) ( 1  -  c~X)  +  2h(A) (1  -  e~X) 


J  h(u)  e  U  du 


+  2h (6)  e"6  h (A) (1  -  e"X)  \  < 

2  -A 

as  long  as  h  ( A )  ( 1  -  e  )  Is  bounded  as  A  -*  0  .  But  this,  of  course,  holds  if 


e 

J " h2(u)  e  U 


du  <  a  ,  and  hi  in  (0,  6)  . 


One  can  dispose  of  the  other  pieces  of  the  integral  by  similar  arguments. 
From  (7.11)  and  (7.12)  it  readily  follows  that 


(s)  -*  a(s) 

n 


(7.17) 


for  every  0  <  s  <  1  and  again  by  (7.12),  (7.16)  and  the  dominated  convergence 
theorem  the  lemma  follows.  Suppose  we  can  show 


lim  sup  Var  R  .  *  0 

M)  -  nl 


(7.18) 


for  i  *  1,  2  .  We  claim  the  theorem  follows  under  our  assumptions.  To  see  this 
note  that  (7.9)  and  (7.18)  imply  that  for  every  e  >  0  ,  there  exists  a  6 
such  that 


lim  sup  d(T  ,  +  yM  < 

rn  \  n  n  n/  - 


(7.19) 


where 


d(X,  Y)  -=  E  { |  X  -  Y  |  (1  +  |  X  -  Y  |)-1} 


(7.20) 


is  the  usual  metric  for  convergence  is  probability.  Now,  (7.3)  yields  for  6  as 
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above , 


lim  sup 


n  d(T  n-i5  l  a ,  (i (n  +  1)_1)(D.  -  1)  4  n^  p6  4  /  ]  < 
n  V  n  6  l  n  n  J  - 


(7.2)) 


and  by  Lemma  7.1,  for  6  sufficiently  small. 


lim  sup  d  |T  ,  S  4  n*5  4  < 

n  \  n  n  n  n  /  - 


(7.22) 


This,  of  course,  implies  that  there  exists  a  sequence  of  constants  such 


that , 


lim  d(T  ,  S  4  K  )  -  0  .  (7.23) 

n  n  n  n 

Since  and  are  both  asymptotically  normal  with  mean  0  by  Lemma  4.1 

of  [2]  and  (ii)  and  the  central  limit  theorem,  it  foilows  that  Kr  -*  0  .  We 
prove  (7.18)  for  i  *=  1;  i  ■=  2  is  proved  analogously.  Let  6  be  defined  as 
in  the  discussion  preceding  (7.15).  Define, 


-  h(Xi)  if  0  <  <  -log(l  -  6) 

■  otherwise 


(7.24) 


where  the  are  uniformly  distributed  on  (-log(l  -  6),  B)  and  independent 

of  each  other  and  the  .  B  is  so  chosen  that  the  density  of  at  h(-log(l 

is  the  same  as  the  density  of  h(X^)  at  rvlog(l  -  5))  . 

We  argue  as  in  (2].  Let  Z.,.  <  ...  <  Z.  .  be  the  order  statistics  of  the 
°  (1)  (n) 


Then, 
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*7  ]  <h(x(i))  -  z<i)> 

\i<6n 


Efi<L (h(X(i)  ~ Z(i))  i -  «>> 


(7.25) 


-2lLv  1 


i <6n  '*>  IX(1)  >  -]og(l  -  6)] 


2  I 


i<6n 
2 


jh  (X(0)  1(X(1)  >  -log(l  -  6)  ]  j 


The  first  identity  follows  by  definition  of  the  ,  the  second  inequality 

from  the  and  Minkowski  inequalities.  Now, 


E/i<Lz«)  ''"a) ; 

2 


2  2 

2  2  6  n  .  i/i 

6))\  -  4  P|X("6n1  >  "l0E(1 


-  6)]  -  0  (7.26) 


since  by  [2]  Lemma  2.2,  P 
On  the  other  hand, 


Xr.  “i  >  -log (1  -  6)  ■*  0  exponentially, 

on  - 


E/h  (X,,J  1 


(i)  lx(i)  >  -logO  ~  «)] 


/  h 2 ( v )  n!  ( (i  -  1)!  (n  -  i) I ]_1  e-(n"i+1)v 


V>-l0g(l-d) 


(7.27) 


(1  -  e~v) 1-1  dv 


- (n-i)v  . 

sup  e  (1 

v>-log(l-<5) 


"V'1  E(h2(X1)) 


Now, 
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Finally, 


sup  e 
v>-log(l-<5) 


- (n-i)v 


“Vvi-1  /n-l\ 

a  -  f 


fi-l  .,n-i  /n-1  \  6 

6  (1  -  o)  n\i-l/  ^or  i  <  ~2  n  • 


(7.28) 


sup  61  1  (1  -  6)n  1  n I*} 
i<  6n  1-1 


[f]-= 


Jlnl 

n— 


■(£) 


(7.29) 


by  an  easy  induction  argument  on  i  .  By  Lemma  2.2  of  [2]  the  right-hand  side 

of  (7.29)  0  exponentially.  Then,  (7.27),  (7.28)  and  (7.29)  imply  that 

^  1 

lim  lim  sup  Var  R  .  ■=  0  if  and  only  if 
6^0  n  nl 


lim  lira  sup  Var|  ]  Z,.j  B  0 

n  Li;®’"  J 


(7.30) 


where  are  defined  for  fixed  6  .  Since  the  Z^'s  are  independently  and 

indentically distributed  with  a  density  positive  an  its  convex  support  and 

E(Z^)  <  «°  ,  (7.30)  follows  from  [2],  The  sufficiency  of  (i),  (ii)  and  (iii)  b) 

for  the  conclusion  of  the  theorem  follows.  The  sufficiency  of  (i),  (ii),  and 

(iii)  a)  is  an  easy  consequence  of  theorem  of  [3], 

6  2 

Using  the  methods  of  [2]  on  T  one  can  show  that  E(S  -  T  )  -*•  0  . 

b  n  n  n 

As  in  [3]  the  smoothness  conditions  on  h  may  be  relaxed. 
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